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Abstract. In the Landau theory of phase transitions one considers an effective 
potential $ whose symmetry group G and degree d depend on the system under 
consideration; generally speaking, $ is the most general G-invariant polynomial 
of degree d. When such a $ turns out to be too complicate for a direct analysis, 
\^ . it is essential to be able to drop unessential terms, i.e. to apply a simplifying 

f^) ' criterion. Criteria based on singularity theory exist and have a rigorous foun- 

dation, but are often very difficult to apply in practice. Here we consider a 
simplifying criterion (as stated by Gufan) and rigorously justify it on the ba- 
sis of classical Lie-Poincare theory as far as one deals with fixed values of the 
C^H' control parameter(s) in the Landau potential; when one considers a range of 

(— I values, in particular near a phase transition, the criterion has to be accordingly 

partially modified, as we discuss. We consider some specific cases of group G as 
examples, and study in detail the application to the Sergienko-Gufan-Urazhdin 
model for highly piezoelectric perovskites. 
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Introduction. 

Landau theory is a powerful and widely used tool in analyzing second order phase 
transitions. This theory was introduced essentially on a phenomcnological basis 
[1, 2], but nowadays it can, in principles, be given a very rigorous justification in 
the framework of singularity theory [3] , or more precisely in that of catastrophe 
(or, in the russian parlance, perestroijka) theory [4, 5]. 

In the same framework, one can also produce constructive techniques to 
simplify the power expansion describing the effective potential at the transition; 
thus "the art of throwing away inessential terms of Taylor series, retaining the 
higher order but physically important terms" (quoted from [4], sect. 5.5) can be, 
in principles, replaced by rigorous computations based on the so called "spectral 
sequence" [3, 7]. 

We have written twice, above, "in principles": indeed, the singularity theory 
approach to Landau theory is rigorous but extremely difficult to implement 
(beside being not so well known among the physicists who apply Landau theory 
to concrete problems), and would require quite involved computations even in 
the simplest case. 

Thus, in practice, what happens is that once the relevant symmetry group 
G has been identified (this is a physical - and not mathematical - matter, 
which we assume to be known by readers), and the general G-symmetric Taylor 
expansion for the effective potential $ has been written down (in terms of G- 
invariant polynomials), one resorts to semi-empirical criteria to simplify the 
expression of $. 

The validity of these criteria is essentially proved a posteriori by checking 
that the simplified potential correctly describes the experimental results. 

These simplifying criteria concern essentially two matters: (a) the order N 
at which the expansion can be truncated; (b) the terms of order smaller than or 
equal to N which are "inessential" , i.e. that can be dropped without changing 
the qualitative properties of the potential $. 

As for point (a), the essential criterion is that of thermodynamical stability: 
we require that there is some neighbourhood T> of the origin which is invariant 
under the evolution determined by $, i.e. such that — V$ points inward at &D. 
In practice, we require that the higher order terms have positive coefficients, so 
that at sufficiently large \x\ the potential $(a;) is convex (see below for definition 
of the notation and a more precise statement). 

In the present paper we want to discuss point (b). A popular criterion for 
this [8] is the following, which we quote almost verbatim from [9] (the term 
"minimal integrity basis" appearing here will be explained below): 

Simplifying criterion. Let G be a compact Lie group, acting in R™ through 
a linear representation; let { J\, ..., J r } be a minimal integrity basis for G, and 
let F(Ji, ..., J. r ) : R" — ► R be a potential. Define, for i = 1 , ..., r and with (., .) 
the scalar product in R™, the quantities 

r dF 

Ui(Jl,...,J r ) := J2 7)T ( VJ fe' VJ • (° J ) 

fc=l k 



A term can be omitted in the potential function without violating the type of 
extremal behaviour of this function if its coefficient is small enough and if it can 
be expressed as 

r 

J2 Qi(Ji,-,Jr) Ui(J lt ...,J r ) + h.o.t. (0.2) 

1=1 

with Qi polynomials in Ji,..., J r and "h.o.t." denoting higher order terms. 

In the present note we want to clarify the meaning of this criterion (see 
also [9, 10]), which should be seen in the context of the orbit space approach 
to variational problems [17, 18, 19, 21, 22], and discuss it in the light of the 
theory of (symmetric) Poincare-Birkhoff normal forms [11, 12, 13, 14, 15]. This 
will enable us to understand in simple terms the origin of the criterion, and 
how it should be modified when considering a range of values for the control 
parameter (s) 

Thus, the goal of the paper is twofold: (i) on the one hand, provide a 
rigorous proof of a reduction criterion (see sect. 3) closely related to the Gufan 
simplifying criterion quoted above when we consider the Landau potential at a 
given value of the control parameter(s); (ii) on the other hand, show that when 
discussing a full range of values for the control parameter(s), we need some 
further restriction. 

It should be stressed that Sanders [16] has recently clarified the relation be- 
tween normal forms, spectral sequence theory, and cohomological issues; how- 
ever, the computations required in this framework are - as for spectral sequences 
in general - extremely involved already in the simplest case. We will stay at a 
simpler although less general level. 

The plan of the paper is as follows. In section 1 we set our general notation 
and recall some general group-theoretical results to be used later on. Section 
2 is devoted to discussing the so-called Poincare, and Lie-Poincare, changes of 
coordinates. In section 3 we apply these to a G-invariant potential; this will 
show that indeed a number of terms - those satisfying a clear-cut condition - 
can be eliminated by a well defined algorithm. The discussion of sections 1-3 is 
conducted at a given value of the control parameter (which thus play little role 
in it) and provides mathematical justification for the reduction criterion, closely 
related to the simplifying criterion given above, stated at the end of section 3. 
In section 4 we discuss how this result applies to Landau theory, taking into 
due account the role of external parameters (such as temperature) controlling 
the phase transition. In section 5 we apply our discussion to a specific case, i.e. 
the Gufan-Sergienko model for highly piezoelectric perovskitcs [9]; we will only 
discuss the form of the Landau potential to be considered, and not the physical 
consequences of it. 
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1 Basic notations and Landau theory 

In this section we assume the reader has some general knowledge of Landau 
theory, as provided e.g. by [2] (see e.g. [6] for more detail), and - as far as this 
is concerned - just set the notation to be used below. 

On the other hand, we need to recall some results concerning G-invariant 
polynomials, with G a Lie group acting (linearly) in a finite dimensional space - 
physically, the order parameter space; a huge mathematical literature is devoted 
to this subject. We will try to keep to as simple an exposition as possible; 
there are expositions designed for physicists rather than for mathematicians 
[17, 18, 19, 21, 22], and the reader desiring more details is referred to these, 
see in particular the first part of [17] and the first chapter of [22]. A readable 
introduction to the mathematical point of view is provided by [31]. 

We denote by x € M ~ R m the order parameter, and by G the group acting 
in M to describe the symmetry of the system under study, or more precisely 
the symmetry inherited in the order parameters space. 

We stress that G acts through a (real) representation, i.e. a set of matrices 
{T gi g € G}; however this is fixed once and for all, so that to avoid cumbersome 
notation we identify g and T g . 

Due to a well known theorem of Palais and Mostow [17], we can always 
assume that the G action is orthogonal after a suitable mapping to a higher 
dimensional space; however, the group actions met in applications of Landau 
theory are generally already orthogonal. We will thus assume that T g are or- 
thogonal matrices; this implies that one of the G-invariants is always |x| 2 . 

The effective potential (fr(x) £ R is a G-invariant polynomial, so that we 
should determine the most general polynomial in x , ...,x m which is invariant 
under G. The Landau potential Q(x) will be a truncation of $ to a suitable 
order TV; moreover, we will be able to omit some "unessential" terms of order 
n < N (see below). 

The polynomials $ and $ have coefficients depending on external parame- 
ters A (e.g., temperature, applied magnetic field,...); the state of the system is 
described by the minima of Q\(x), which we denote as x a (X). Note that, in 
general, there will be different minima for a given value of A; in particular, if 
x(A) is not a fixed point for the G-action, then the whole G orbit through x(X) 
will be an orbit of minima (if G acts effectively in M, then by definition is the 
only fixed point for G) . Moreover, there can be different G-orbits of minima for 
a given value of A. 

The symmetry of the state corresponding to x(X) will correspond to G x /\\, 
the isotropy group of a; (A). We recall that by definition 

G x := {g e G : T g x = x } . 



Note that if x and y are on the same G-orbit, then y = T g x for some g € G, and 
then G y = gG x g~ x : there is a conjugacy class of isotropy subgroups associated 
to any G-orbit in M, and this is also called the orbit type [Gx\. A phase will be 
described precisely by an orbit type. 

Example 1. Consider G = Z2 = {e,g} acting in R via gx = —x. Then 

Go = G, while for any x ^ 0, Ga; = {e}. An invariant Landau potential is of 
the form $(x) = — Xx 2 /2 + x 4 /4 (with A <G R), and as A changes sign we pass 
from a symmetric phase (orbit type G) to a non symmetric phase (orbit type 

M). " 

Example 2. Consider G = SO(n) acting in R™ through the defining repre- 
sentation. Then any point x ^ is left fixed by SO(n — 1) acting as rotations 
around the axis identified by the origin and x. All these subgroups are mapped 
one into the other by rotations in SO(n). When we consider the potential 
$(2;) = — Ax 2 /2 + x 4 /4 (A <G R) we observe a phase transition as A changes sign, 
and in this language it corresponds to passing from a phase with orbit type 
[G{0}] = SO(n) to one of orbit type [Gx] = SO(n - 1). 

In order to have a phase transition at a value A = Ao, it is needed that the 
orbit type [Gx(A)] is not constant in a neighbourhood of Ao, no matter how 
small. 

In discussing what the "suitable" order N is and which terms of order n < 
N can be omitted, we need to introduce some notions relative to G-invariant 
polynomials. 

1.1 G-invariant polynomials 

We will consider a compact Lie group G acting linearly and orthogonally in the 
space R" (many of the notions and results mentioned below have a much wider 
range of applicability, but here it suffices to consider this frame on the basis 
of physical needs). We look at the ring of G-invariant scalar polynomials in 
x 1 ,...,^™, denoted as S(G). 

By the Hilbert basis theorem [17, 24, 25], there is a set {J\(x), ..., J r (x)} 
of G-invariant homogeneous polynomials of degrees {d\, ...,d r } such that any 
G-invariant polynomial $(x) can be written as a polynomial in the {Ji, ..., J r }, 
i.e. 

$(x) = ^[Jx(x),...,J r (x)] (1.1) 

with 7? a polynomial in (Ji , ..., Jr). 1 

With our hypotheses on G, the algebra of G-invariant polynomials is finitely 
generated, i.e. we can choose r finite. When the J a are chosen so that none 
of them can be written as a polynomial of the others 2 and r has the smallest 
possible value (this value depends on G), we say that they are a minimal integrity 

J It may be interesting to note that the theorem also holds for smooth functions [28]. 
2 It is essential to understand that some of the J a could be written as non-polynomial 
functions of the others, and the J a could verify polynomial relations, sec Example 5 below. 



basis (MIB). In this case we say that the { J a } are a set of basic invariants for 
G. There is obviously some arbitrarity in the choice of the h in a MIB, but the 
degrees {di, ..., d r } of { Ji, ..., J r } are fixed by G. 3 

We will from now on assume we have chosen a MIB, with elements { Ji, ..., J r } 
of degrees {d±, ..., d r } in x, say with di < d<i < ... < d. r . 

Example 3. The simplest example is that of the fundamental representation 
of orthogonal groups themselves, G = 0(m), acting in M — R m . In this case 
(recall that invariants separate orbits) the only invariant function is r — |x|; 
however as we have to deal with polynomials in the coordinates (x , ..., x n ), the 
basic invariant will be J\ = r 2 . ^> 

Example 4. Consider R 2 with coordinates (x,y), and in this the group G 
generated by the elements g x , g y acting as g x : (x, y) — > (—x, y) and g y : (x, y) —* 
(x, —y). Basic invariants for this are obviously J\ — x 2 and J2 = y 2 ■ (} 

Example 5. Let us consider another simple example: M — R 2 , and G is the 
group generated by simultaneous reflections in x and in y, i.e. G = {/, (R x R y )} 
with R x : (x,y) — > (— x,y) and R y : (x,y) — » (x, — y). It is immediate to sec 
that any G-invariant polynomial can be written as $(x,j/) = 7?(x 2 ,y 2 ,xy), i.e. 
we have a basis of three quadratic polynomials: J\(x,y) = x 2 , J2(x,y) = y 2 , 
and Js(x, y) = a;y. Note that here no elements of the MIB can be written as an 
algebraic (that is, polynomial) function of the others, but they are nevertheless 
satisfying an algebraic relation: indeed, J1J2 = J 2 . That is, they are not 
algebraically independent. <) 

Example 6. Let us consider another simple example (relevant to the situation 
considered in [9]): now M — R 3 , with G the group generated by independent 
reflections in x, in y and in z (so, with the notation used in Example 5, G 
is generated by R x , R y and R z ). We have now a basis of three polynomials 
J a (x,y, z), given by J\ = x 2 + y 2 + z 2 , J 2 = x 2 y 2 + y 2 z 2 + z 2 x 2 , J 3 = x 2 y 2 z 2 . 



When the elements of a MIB for G are algebraically independent, we say 
that the MIB is regular, if G admits a regular MIB we say that G is coregular. 
This case is, needless to say, easier to analyse; as the Example 5 above shows, 
the non-coregular cases can include very innocent-looking groups. 

An algebraic relation between elements J a of the MIB is said to be a rela- 
tion (or, with a term of astronomical origin, a sygyzy [23]) of the first kind. The 
algebraic relations among the J are a set Sru of polynomials in the { J\, ..., J r } 
which are identically zero when seen as polynomials in x; this admits a a min- 
imal set of homogeneous generators {s[ , s 2 , ..., sL'n}. If there are algebraic 
relations among these, they are called relations (or syzygies) of the second kind, 
and so on. A theorem by Hilbert guarantees that the chain of syzygies has 

3 They are determined through the Poincare series of the graded algebra Pq of G-invariant 
polynomials, sec [19]. 



finite maximal length (this is the homological dimension of the graded algebra 
Pg mentioned above). 

In the following we will need to consider a matrix built with the gradients 
of basic invariants, which we call (following Sartori [18, 19]) the P-matrix. This 
is defined as 

V ih {x) := (VJi{x),VJ h {x)) (1.2) 

with (., .) the standard scalar product in M = R m . 

Note that the gradient of an invariant is necessarily a covariant quantity; the 
scalar product of two covariant quantities is an invariant one 4 , and thus can be 
expressed again in terms of the basic invariants. Thus, the V -matrix can always 
be written in terms of the J themselves. 

Example 4 (continued). The basic invariants J\ — x 2 and J 2 = y 2 are inde- 
pendent. Their gradients are VJi = (2x, 0), \7J 2 = (0, 2y); the P-matrix is 
hence 

v = (4x 2 \ /4Jx 

V Ay 2 ) \ 4J 2 





Example 5 (continued) . If J\ — x 2 > 0, J2 = y 2 > and J3 = xy, then as 
(x, y) varies in R 2 , the point (Ji, J2, J3) € R+ X R+ x R c R 3 varies on the 
manifold J 2 = J\ J 2 - 

The gradients of the basic invariants are VJi = (2a;, 0), VJ2 = (0,2j/), 
^^3 = {Vt x )\ note there is an obvious linear dependence among these at all 
points, which shows the J arc functionally dependent. The "P-matrix is 



4x 2 





2xy \ 


/4Ji 





2J 3 





4y 2 


2xy 


= ° 


4J 2 


2J 3 


2xy 


2xy 


x 2 + y 2 I 


\2J3 


2J 3 


J1 + J2 



P(x,y) 



Example 6 (continued). In this case we have 

/(j/ 2 + z 2 ).t\ /(y 2 z 2 )x* 

, VJ 2 = 2 (a: 2 + z 2 )y , VJ 4 = 2 (a; 2 z 2 )2; 
\(x 2 + y 2 )z) \(x 2 y 2 )z, 

Therefore, with straightforward explicit computations, the P-matrix is 

2J 2 3J 3 

P = 4 ( 2J 2 (J1J2 + 3J3) 2J!J 3 
2J1J3 J2J3 



♦ 




♦ 



4 Had we been considering a general manifold M with metric gij, the scalar product be- 
tween covariants quantities would be defined through g 1 ^ , i.e. we would have 7'i/ l (x) = 
(yj i ) a g<*f>(VJ h ) l 3. 



1.2 Application to Landau theory 

Let us briefly comment on how the mathematical results mentioned above are 
relevant in Landau theory. As implied by the Hilbert basis theorem, any G- 
invariant polynomial, such as the effective potential $(x) or the Landau poten- 
tial ^(x), can be expressed as a polynomial tt(J(x)). 

In this way the evaluation of the map $ : M — >R is in principles substituted 
by evaluation of two maps, h : M — >Q and n : il — >R; here wc have denoted 
by il C R r the target space for J = (Ji, ..., J r ). 

However, if - as is the case in Landau theory - we have to consider the most 
general polynomial in S(G), we only have to deal with the map 7r : Q — > R. 

Note that for corcgular groups we have simply fi = R r , while for non- 
corcgular groups Q is a submanifold of R r . 

The space il is also known as the orbit space for the G action on M; indeed, 
its points are in one-to-one correspondence with the G-orbits in M. 

In facts ft is a semialgebraic manifold, i.e. a submanifold in R r defined by 
algebraic equalities and inequalities; it is moreover a stratified manifold, i.e. 
the disjoint union of smooth manifolds of possibly different dimensions, with 
bordering relations related to the conjugacy class of isotropy subgroups on the 
G-orbits represented by points on each of these manifolds. We will not discuss 
these points, but just refer to [17, 18, 19, 21, 26]. 

As the J are nontrivial homogeneous polynomials, the origin of R m will cor- 
respond to Ji = ... = J r = 0, and a neighbourhood of zero in M will be mapped 
to a neighbourhood of zero in fi; thus we will not have to actually evaluate the 
map J : M — > fl, but can just focus on 7? : ft — > R in a neighbourhood of the 
origin. 

We will boldly summarize our discussion as the following 5 

Landau-Michel principle. 

Landau theory can be worked out in the G-orbit space tt :— M/G. 

1.3 Thermodynamic stability 

At this point we can already briefly discuss how the request of thermodynamic 
stability, i.e. convexity (see the introduction) is reflected in the polynomial 

Let us first of all consider the coregular case; now 7? : R r — > R, and the J a 
can be considered as independent variables. The minimal Landau polynomial 
$(x) = 7r( J) will be quadratic in the J, and the stability is ensured by requiring 
that the matrix of second derivatives Dj/, = d 2 n/(dJidJh) is positive definite. 

So the prescription in this case will be to consider a polynomial of order 
N = 2max((ii, ..., d r ) = 2d r (recall d\ < c?2 < ••• < d r ); and of course choose 
coefficients so that the matrix D is positive definite for large |x|. 

5 Louis Michel (1923-1999) pioneered the use of orbit space techniques in Physics and 
Nonlinear Dynamics, originally motivated by the study of hadronic interactions; see e.g. [22, 
26, 27]. For an essential bibliography of his works on this theme, see [38]. 



Note that if we deal with a non-corcgular case this prescription also works: 
maybe it would be possible to stop at a lower order, as we have to care only 
about the submanifold of fi allowed by the relations between the J a , but if we 
require stability in all of fi we are on the safe side. 

We stress that the prescription is not to write 7r as a quadratic polynomial 
in the J a and then express $ in terms of this; rather it is to consider the most 
general G '-invariant polynomial of order 2d r . This can contain quite high powers 
in some of the J a 's, see example 6. 

It is essential to recall that the coefficients of (at least some of) the poly- 
nomials will depend on the external parameters; in particular, this will be the 
case for J\ — |x| 2 , whose coefficient controls the loss of stability of the critical 
point x = and thus the onset of the phase transition. 

Example 4 (continued) . In the case of example 4, the basic invariants are both 
quadratic, so we consider an expansion quadratic in these. That is, we consider 

$ = $i + $ 2 ; 

<I>i = aiJi + a 2 J2 , 

*2 = (biJi+b 2 Ji) + ciJiJ 2 ■ 

Thus <& depends on 5 coefficients. The origin is stable provided both the a% are 
positive; thus (at least one of) the a& should depend on parameters, and pass 
through zero, in order to have a phase transition. This example will be studied 
in section 4 below. ^> 

Example 5 (continued) . In the case of Example 5, all basis polynomials are 
quadratic in (x,y), so we have to consider again an expansion up to order four 
in (x,y), hence quadratic in the J a . This means we will consider 



$ = 


= $1 + $2 5 


$1 = 


= {a\J\ + a 2 J 2 + C13J3) , 


$ 2 = 


= [(biJl + b 2 Jl + 63 4) + (c^Js 



C 2 JlJ 2 )\ ■ 

We have not included terms of the form cqJiJ 2 in $2: as J\ J 2 = Jf, these can 
be absorbed in &3</f . Hence $ depends on 8 coefficients; note that again (at 
least some of) the at should depend on parameters to have a phase transition. 
This example will also be studied in section 4 below. ^> 

Example 6 (continued). In the case of Example 6, we have d\ — 2, d 2 = 4, 
(I3 — 6. Hence we have to consider an expansion up to order 12 in (x, y, z). We 

write 

$ = $1 + $ 2 + $3 + $4 + $5 + *6 

and the homogeneous terms $ m (of order 2m in x, y, z) are given by 



$1 


= a\3\ 


$ 2 


= a 2 J 2 + biJl 


$3 


= a 3 J 3 + b 2 jf + ci Ji J 2 


$4 


= b 3 Jf + b 4 J% + C2J1 J 2 + C3J1J3 


$5 


= 64 Jf + C 4 jfj 2 + C 5 jfj 3 + CqJiJ^ + c 7 J 2 J 3 


$ 6 


= Mi + K J 2 + b 7 Jl + c 8 JfJ 2 + c 9 jfj 3 + CioJfJt 



CnJiJ 2 J3 



Thus <& depends on 21 coefficients; it is obvious that we need to simplify it in 
order to be able to discuss its behaviour. 

In [9], Sergienko, Gufan and Urazhdin argued that in this case one can just 
consider 

$ = aiJi + (3i J 2 + 71^3 + a 2 Ji + /3 2 J| + 72 -/f , 

i.e. a potential quadratic in the basic invariants and with no mixed terms (terms 
of the type JiJ k with i ^ k); this was further justified and considered in [9]. 

Note that the stability of the origin is controlled by the sign of the coefficient 
a\. so this should depend on parameters, and pass through zero, in order to 
have a phase transition. 

This example will be studied in detail in section 5 below. <0> 



2 Poincare and Lie-Poincare transformations 

In this section we recall a technique which is the fundamental tool of the theory 
of Poincarc-Birkhoff normal forms [11, 12, 13, 14, 15], i.e. Poincare transforma- 
tions and an improved version of these, Lie-Poincare transformations. 

Were we dealing with the complete Taylor expansion $(x), we would need 
the whole normal form theory; as we only consider a truncated polynomial $>(x), 
we will only need some part of the theory. 

Consider a neighbourhood D of the origin in R" with coordinates (x 1 , ..., x n ); 
and the Taylor series for a smooth function F(x) G R, vanishing in the origin: 

oo 

F(x) = ]TF fc (x) ; F k (ax) = a k+1 F(x) . (2.1) 

fc=0 

Note that F k is homogeneous of degree k + 1: the reason for this notation 
will be clear below. We will denote the set of smooth scalar (vector) functions 
homogeneous of degree k + 1 as S k (14). 

We want to consider near-identity changes of coordinates in D generated by 
a homogeneous vector function h k <E V k with components h l m £ Sk (Poincare 
transformations), and how these are reflected in the expansion (2.1) for F. 

If we write the change of coordinates in the form 

x l = y l + K m (y) ■ hl n eS m (2.2) 

then we have (using an obvious simplified notation) 

n f)F n 1 r) 2 F 

F k (x) = F k (y + h)= F k (y) + £ *tf + £ - g-^hW + ... (2.3) 

j=l y i,j=l y y 

Here the term with with q derivatives (q — 0,1,...) belongs to S k + qm , with 
9 = 0,1, 
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This shows at once how the F k are transformed; if we write F(x) = F(y) 
"Ylik F k (y), then under (2.2) we have 



F k -- 
F k -- 


= F k 
= F k 


+ (h m 


■V)F fc 


F k -- 


= F k 


+ (h m 


•V)F fc 



for k < m 
for m < k < 2m 
+ {l/2)h m AF k _ 2m for 2m < fc < 3m 



(2.4) 



A modification of this technique, more suited to our needs, is known as Lie- 
Poincare (LP) transformations [32, 33, 13, 35]. With this, rather than consider- 
ing the change of coordinates (2.2), we will consider the change of coordinates 
obtained as the time one flow of the vector field x % = h l m (x): if $(i; Xo) is the 
flow under this with initial condition x(0) = xq then we define the change of 
coordinates as x = $(l;y), also written as x = e y. 

Needless to say, this just coincides with (2.2) at first order; however, this 
approach has several advantages. Leaving apart the theoretical ones [35], we 
stress that LP transformations are better suited to deal with the case where 
the function F has some symmetry property which we want to preserve in the 
change of coordinates. 

The last statement is better understood, anticipating slightly our discussion, 
considering the gradient / = VF; we write f = ^2 fk, with f k = WF k+ i <G V k . 

Under a LP change of coordinates we have (see e.g. [14, 29] for details) 

[fe/m] 
f" = E - nS m{fk-sm) ; (2.5) 

here we have denoted by [a] the integer part of a, and introduced the operator 
H m defined by 

n m (f) := {h m , /} := (hm ■ V)/ - (/ • V)h m ■ (2.6) 

Note that (2.5) is just the Baker-Campbcll-Haussdorf formula 6 ; it states that 
the vector function f(x) is written as f(y) with 

/ = [e th fe~ th ] t=1 . (2.7) 

Let us now briefly comment on how these changes of coordinates are used 
(we keep to the vector case for ease of discussion); this rests on (2.4). 

With both the Poincare and Lie-Poincare approaches, if we consider a change 
of coordinates with generator h m € S m , then: (a) the terms of degree k < m are 
not changed at all; and (b) the terms of degree m (actually, of degree m < k < 
2m.) are changed in a very simple way, depending only on the linear part of the 
system, i.e. as F m — F m +H m (fo). (Higher degree terms change in a slightly 
different way depending if we are using the Poincare or Lie-Poincare approach) . 

6 Here the role of the commutator is taken by the bracket {.,.}, which is nothing else than 
a translation of [., .] in terms of components of the commutator \Xh,Xf\, with X^ = h l di and 
Xf = f l di. Indeed, for X^ and Xf as above, [Xh,Xf] = {h, f} l di. 
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Property (a) allows to consider sequentially changes of coordinates with gen- 
crating functions homogeneous of degree 2, 3, ....; each time one is not changing 
terms of lower degree - so that terms in Sk can be considered as stable after 
step k - and by property (6) one is able to eliminate all terms in the range 
of the operator £ : = {/o, •}• Terms in the complementary space cannot be 
eliminated, and are said to be resonant. 

Hence, by applying repeatedly changes of coordinates of this type, one can 
reach a system which is in normal form - namely, contains only resonant terms 
- up to any desired order N. 

3 Transformation of invariant polynomials 

3.1 Lie-Poincare transformations and invariant polynomi- 
als 

Let us apply the technique illustrated in the previous section to G-invariant 
polynomials <&(x) = ir[J(x)]. We write 

oo 
fe=0 

where ^k(ax) = a k+1 ®(x). 

We want to consider changes of coordinates of the form 

.t s = y* + h\y) ; (3.2) 

moreover, we want to preserve the symmetry properties of $: the change of 
coordinates should be G-equi variant. Thus the function h : M — ► M has to 
transform in the same way as y under the G-action, i.e. we have to require that 

h(T g y) = T g h{y) \/y € M \/g € G . (3.3) 

We will choose h to be the gradient of a G-invariant function H(x), i.e. 7 

h\y) = 5- d JM . (3.4) 

We stress that in general there can be equivariant vector polynomials which 
are not obtained as the gradient of invariant scalar polynomials 8 ; however, the 
choice (3.4) guarantees that (3.3) is satisfied: we are not as general as possible, 
but we are on the safe side. 



7 For general M with metric g, the 8 in the formula (3.4) should be replaced by g. 

8 This can be the case for very simple groups: it happens e.g. for any representation of 
SO (3) at the exception of the fundamental one, and h of degree higher than 3; see [30] for 
details. 
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In order to know how (3.2) acts on (3.1), it suffices to know how it acts on 
the basic invariants J a . This action is simply 

J a (x) = J a [y + h{y)] = J a {y) + {5J a ){y) 

= J a (y) + (d,J a /dyP)hP + \(d 2 J a /dyPdyi)hPW + ... ^° J 

and the action on $ is therefore 

$(x) = w[Ji{x),...,J r (x)] = 

= 7r[Ji(y),...,J r (y)]+EU(*r/^a)[(*Ja)(y)] (3.6) 

+ (1/2) E^i (d 2 n/dJ a dJ b ) [(SJ a )(y)] [(5J b )(y)] + 

Note that if $ is of finite order TV, this change of coordinates will in general 
produce terms of higher order, i.e. $ will be changed into a polynomial of order 
higher than N. Thus, the truncation to order TV = 2d r should be performed 
again after all the required Poincare changes of coordinates took place (or after 
each one, for computational convenience). 

The formulas (3.5), (3.6) can be considerably involved, and in principles can 
be read off equation (2.5). Luckily, for our discussion we only need the first 
order terms; dropping higher order terms, (3.5) reads 

J a (x) = J a {y) + (dJ a /dy p )hP ; (3.7) 

recalling (3.4), this is 

J a (x) = J a (y) + (dJ a /dyF)5 pq (dH/dyi) . (3.8) 

As H is G-invariant, it is also possible to write it as a function of the basic 
invariants: H{x) = x[Ji(x), ..., J r (x)]. Hence, 



dH_ _ dH_ <9Jb 

dy k dJb dy k 



and (3.8) reads 



J a (x) = J a (y) + (dJ a /dyP)5P«(dH/dJ b )(dJ b /dyi) 

= J a {y) + ELi(v-/a|vj 6 )09tf/a/ fc ) . 



(3.9) 



(3.10) 



If now we recall (1.2), the above can be rewritten as 

Ja(x) = J a (y) + (SJ a )(y) \6J a := V a fj(dH/dJ b ) . (3.11) 

3.2 Elimination of terms in invariant polynomials 

Let us now apply the above discussion to the reduction of an invariant polyno- 
mial $(x) = ^(Ji, •••, J r )- We have in general 

*(J) -» *(J + 6J) - *(J) + V ^£D-SJ a + h.o.t. . (3.12) 



a=l 
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Under a change of coordinates of the form (3.2), with h given by (3.4), the 
J a (x) change according to (3.11). We will write D a := d/dJ a , and understand 
summation over repeated indices is implied; disregarding higher order terms and 
using the explicit form of (3.11), we get from (3.12) that 

a* dH 

s * = 7rr v ^irr = (D a ^)p aP (D p H) . (3.13) 

OJ a (J<J(i 

Let us now consider the expansion (3.1) for $, and hence for \& = ~^2 k ty k , 
where <$> k (x) := * fe (J). 

Here we are introducing a grading with respect to the x coordinates; note 
that J a is of degree d a , hence V J a is of degree (d a — 1), and the element V a p of 
the "P-matrix is of degree (d a + dp — 2). We denote by dh = m + 1 the degree 
of the function h (which is dH — 1, with dH = m + 2 the degree of the scalar 
function H m ). Thus the element (D a ^ k )V a /3 (DpHm) (no sum on a,/3) is of 
degree 

(k + l-d a ) + (d a +d f3 -2) + (d H -d fj ) = k + d H - 1 = k + m + 1, 

i.e. (D a y k )T a p(D fj H m ) e S m+k . 

This means that under a change of coordinates (3.2) generated by H m the 
terms ^ k with k < m are not changed, while the terms \& m + p change according 
to 

* m+p -^^ m+p = ^ m+p + {D a ^ p )V afi {DpH m ) + h.o.t. . (3.14) 

We can then proceed as sketched in section 2, i.e. operate sequentially with 
H in S\, S2, ■■■■ ; at each stage (generator H m ) we are not affecting the terms 
"3/fc with k < m. Moreover, we can just consider the first order correction, as 
higher order terms are generic and will be taken care of in subsequent steps. 

Consider now, to fix ideas, the case where d± — di = ... = d s = 2 (recall 
this implies there are no invariants of order less than two). Then \I/ TO +i changes 
according to 

* m+ i -» $ m+ i = * m+ i + (D a *i)^(D^ m ) + h.o.t. . (3.15) 

If Wi 7^ 0, we disregard higher order terms; note that necessarily ^1 = a k J k 
(the sum on k runs from 1 to s). In this way we obtain that all terms which are 
of the form 

s r 
fe=l/3=l 

for H G S m +i can be eliminated by our procedure. Needless to say, if the a k 
are all zero we need to consider higher order terms. 

Note also that in order to perform this procedure, we should determine suit- 
able generating functions H m ; this can be done by requiring as many component 
of \& m +p as possible to vanish, and inverting (D a 'fy p )T'ai3 to solve (3.14) for H m . 
The relevant point is, that this requires D a ^/ p ^ 0, and this condition could fail 
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for certain values of the parameters denning * p in terms of the basis invariant 
polynomials. 

Another relevant point should also be addressed: eqs. (3.13) and (3.16) say 
that, in this context (see the discussion below), the functions Qi appearing in 
(0.2) should be identified with the derivatives D^-ff* for a suitable function H* 
(this will be obtained from the generating functions H m used at each step). 
However, this identificiation leads to a natural "compatibility condition" be- 
tween the QiS: that is, we have to require 

DiQj = D,Q t . (3.17) 

Example 4 (continued) . In the case of example 4, the two basic invariants 
J\ = x 2 and J 2 — y 2 are both quadratic, and thus contribute to the first order 
formulas (3.13), (3.15). 

Let us consider a quartic generating function H 2 = hi J 2 + h 2 J 2 + fci J\ J 2 ; 
the quartic term ^ 2 = b\J 2 + b 2 J 2 + cJ\J 2 will change to ^ 2 + <5*2, and from 
our formulas we get 

(5* 2 = (8ai/ii)Ji + {8a 2 h 2 )J 2 + [4(ai + a 2 )k 1 ] ,hJ 2 • 

Therefore, provided the conditions 

ai 7^ , a 2 7^ , ai + a 2 7^ 

are verified, we can choose h\,h 2 , fci such that ^2 = ^2 + 6^2 = 0. These are 
given explicitly by 

h bl u b2 1 c 

hi = , h 2 = , k\ = . 

8ai ~ 8a 2 4(ai + a 2 ) 

Note that actually we do not want to set the term \&2 to zero, in order to 
satisfy the requirement of thermodynamic stability. Thus this is just an example 
of what could be done by the Poincare procedure, but does not apply directly 
to Landau theory. ()■ 

Example 5 (continued). In the case of Example 5, again all the basic invariants 
are quadratic in (x,y). Thus all of them contribute to the first order formulas 
(3.13), (3.15). Recall in this case the invariants satisfy the relation J 2 = J\J 2 . 
We consider a quartic generating function, which we write as 

H 2 = h\Jl + h 2 3l + h 3 J 3 2 + fci Ji J 3 + k 2 J 2 J 3 . 

Under the change of coordinates generated by this, the quartic term in the 

t2 
*2 = *2+<5*2- 



Landau polynomial, ^> 2 — b\ J 2 + b 2 J 2 + 63 J 2 + c\J\J 3 + c 2 J 2 J 3 , will change to 



It follows from our general formulas, with standard algebra and using the 
elation J 2 = J\J 2 , that 

5*2 = {8aihi + a 3 kx)J 2 + {8a 2 h 2 + a 3 k 2 )J 2 + [4(01 + a 2 )h 3 + 3a 3 (fci + k 2 )]J 2 - 



-[4a 3 /ii + 2a 3 h 3 + (6ai + 2a 2 )fci] Ji J 3 + [4a 3 /i 2 + 2a 3 h 3 + (2ai + 6a 2 )k 2 ] J 2 J 3 
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Thus, we can choose {hi, ki} such that \&2 = 9 if the matrix 



/8ai 


a 3 \ 





8a 2 a 3 





4(ai + a 2 ) 3a 3 3a 3 


4a 3 


2a 3 (6ai + 2a 2 ) 


V 


4a 3 2a 3 (2ai + 6a 2 ) / 



has nonvanishing determinant. In fact, with obvious notation, ^ 2 = M(h, k) + 
(b, c). It results 

\M\ = 256[12afa 2 -3a 2 1 a 3 ! + a 2 af + a?(52a 2 ! - 3a|) 

+af (52a^ - 17a 2 a^) + ai(12a| - lTa^af + af)] . 

Let us discuss the vanishing or otherwise of this. If any two of the a, are 
zero, then \M\ — 0. If 01 = (and of course 02 / / a 3 ), then \M\ = 
{d2a\){a\ — 3a|), i.e. we can still set ^ 2 = provided a 3 7^ ±V3a 2 . Similarly, 
if a 2 = (and of course a\ 7^ 7^ a 3 ), then |M| = (aia 3 )(ag — 3a^), and we can 
set ^2 = provided a 3 7^ ±\/3oi- For a 3 = (and of course a\ 7^ 7^ a 2 ), we 
get 

|M| = 4aia 2 (3aJ + 13a^a 2 + 13aia2 + 3a 2 ) ; 

this is nonzero provided a\ 7^ — a 2 , ai 7^ — 3a 2 , a\ 7^ — a 2 /3. 

Finally, if the a{ are all nonzero, the determinant is nonzero provided all of 
the following conditions are satisfied: 

ax 7^ -a 2 

a| 7^ 4aia 2 

a| 7^ 3aJ + 10aia 2 + 3a| 

We won't pursue the analysis of the degenerate cases; the point here is to show 
that all the computations can be done in complete detail using only standard 
algebra. () 

Example 6 (continued). In the case of Example 6, only 3\ is quadratic, so first 
order formulas involve only this, and cannot be used if a\ = 0. 

We will postpone discussion of this example until sect. 5, where we consider it 
in the more general case of varying control parameters, i.e. varying coefficients 
at, hi, a. 

We stress that in examples 4 and 5, we could use the simplified formulas 
(3.13), (3.15), based on linear action (quadratic generating function). However, 
in general we should consider the full change of coordinates generated by H. 
Note that in this case the transformation eliminating terms of a given order will 
at the same time generate terms of higher orders. 

9 Again, one does not really want this in the context of Landau theory, as it would violate 
the requirement of thermodynamic stability; see the previous example. 
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3.3 The simplifying criterion, and discussion. 

The discussion conducted in this and the previous section provides rigorous 
ground for the symplifying criterion stated in the Introduction and due to Gufan. 

It should be stressed that here we have worked with a given Landau poly- 
nomial, i.e. with fixed values of the parameters entering in it. We stress that 
these parameter - or at least some of them - will in general depend on the 
external "control" parameter (temperature, pression, magnetic held, etc), and 
indeed they have to change with these for a phase transition to take place. 

Thus our discussion so far provides a proof of the Gufan simplifying criterion 
only when we work at a given value of the control parameter(s); we will see in 
the next section that some extra care is needed if we want to work on a full 
interval of values of the control parameter(s). 

This is reflected in our statement of the simplifying criterion, see the Intro- 
duction, as there we consider a given potential with no mention of dependence 
upon external parameters. 

Another relevant remark should be made. In this section, and in the whole 
paper, we consider just changes of coordinates: that is, we eliminate terms by 
choosing suitable coordinates, but we are not changing the potential; on the 
other hand, in singularity theory one allows changes of the potential, provided 
these do not alter its qualitative behaviour [3]. Thus, we are not considering 
the most general transformation of $ allowed by Landau theory. On the other 
hand, the transformations considered here are surely allowed, and actually can 
be easily implemented (algorithmically) via a symbolic manipulation language. 

In view of the above, it is worth stating a weaker form of the simplifying 
criterion given in the introduction. This is the result being actually proven by 
our discussion so far; note that it should be called a "reduction" - rather than 
"simplifying" - criterion: indeed, as already remarked, we do not simplify the 
potential, but just provide a reduced expression of it by using more convenient 
(local) coordinates. 

Reduction criterion. Let G be a compact Lie group, acting in R™ through a 
linear representation; let {J\, ..., J r } be a minimal integrity basis for G, and let 
F(Ji, ..., J r ) : R n — ► R be a potential. Define, for i = 1, ..., r and with (., .) the 
scalar product in R™, the quantities 

r dF 
Ui(J u ...,J r ) := V — (VJ,,VJ 4 ) . (3.18) 

t[ dJk 

Let B C R™ be a sufficiently small neighbourhood of the origin in R". Then 
there is a sequence of Poincare changes of coordinates in B, such that the po- 
tential F is expressed in an equivalent form F, where terms of F which can be 
written as 



a=l 



,(Ji,...,J r ) U a {J u ...,J r ) + h.o.t. , (3.19) 



17 



with Q a polynomials in Ji, ..., J r satisfying the compatibility condition 

dQp _ dQ a 
dJ a dJp 

and "h.o.t." denoting higher order terms, can be eliminated in F. 



(3.20) 



The "sufficiently small" in the statement above should be meant in the sense 
that the overall Poincare change of coordinates described by the combination 
of the different ones at each order k = 2, ..., 2d r should dchnc a convergent se- 
ries. In general - for a given finite order d r - this will be the case only in some 
neighbourhood B of the origin. In this respect, it should be stressed that the 
sutuation considered here is naturally characterized by a symmetry, hence the 
symmetric theory (due to Markhashov, Bruno, Walcher and Cicogna) should 
be used 10 . Note however that when one performs a concrete computation, the 
radius of convergence will always be immediately read off the concrete (denom- 
inators in the) inversion formulas. 

4 Reduction criterion with varying parameters 

In Landau theory, one is considering (the vicinity of) phase transitions; that is, 
the coefficients of the polynomial $(x) depend on external control parameters 
A, and necessarily pass through critical values. 

In this section, we investigate how the discussion given so far should be 
modified if we want to consider not just given fixed values of the parameter(s), 
i.e. of the coefficients appearing in the Landau polynomial, but a full range of 
values, including in particular critical ones. 

Let us, to fix ideas, consider A <E R and let J\ — \x\ 2 be the only quadratic 
basic invariant, so that $ = Ci(A) J\ + C2(\)J\ (with 02(A) > 0) and the loss of 
stability of the x — critical point is controlled by the sign of the coefficient 
Ci(A) of J\ in $ - say Ci(A) = —A/2 with standard normalization - so that the 
transition occurs at A = 0, with the symmetry-breaking phase appearing for 
A > (this is e.g. the case for $ = -A|x| 2 /2 + |x| 4 /4). 

As we want to describe a small but finite interval of values of A e A := 
[e_,e+], with e_ < and e + > 0, we have to require that the near-identity 
changes of variables considered in previous sections are defined uniformly in A. 
This means that they must be well defined also at A — 0; in particular, changes 
of variables requiring a division by ci(A) are not allowed. 

More generally (i.e. for general degrees of polynomials), we are allowed to 
consider only those Poincare and Lie-Poincare transformations which are smooth 
and well defined in a full neighbourhood of the critical point, and in particular 
at the critical point itself. 

This means that the reduction criterion stated in the previous section should 
be suitably restricted, as follows. 

10 We will not discuss this theory here, and just refer e.g. to the discussion given in [36, 37]. 
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Reduction criterion with varying parameters. Let G be a compact Lie 
group, acting in R™ through a linear representation; let {J\, ..., J r } be a minimal 
integrity basis for G, and let F(X; J\, ..., J r ) : A x R™ — > R be a potential 
depending on the control parameters A G A C R p . Define, for i = 1, ...,r and 
with (., .) the scalar product in R", the quantities 

r dF 
t^(A;Ji,...,J r ) := V j— (VJ fc ,VJ 4 ) . (4.1) 

Let B C R™ be a sufficiently small neighbourhood of the origin in R". Then 
there is a sequence of Poincare changes of coordinates in B, such that the poten- 
tial F is expressed in an equivalent form F in a neighbourhood A of the critical 
value A = Ao in the space of the control parameters X; terms of F which can be 
written uniformly in A, as 

r 

Y, Qi(Ji,-,J r ) l7i(A;Ji,...,J r ) + h.o.t. (4.2) 

8=1 

where Q a are polynomials in Ji,...,J r satisfying the compatibility condition 
{dQ a /dJfj) — {dQfj/dJ a ) and "h.o.t." denotes higher order terms, can be elim- 
inated in F . 

Note the only difference with respect to the criterion given in the previ- 
ous section lies in the dependence on A and in requiring uniformity in A, as 
emphasized in the statement. 

In order to better understand this point, it is useful to consider again in 
detail the elementary example 4; we will then also analyze in detail example 5, 
while example 6 is analyzed in the next section. 

Example 4 (continued). In the case of example 4, a quadratic generating 
function H\ would produce a linear change of coordinates (we are not interested 
in these) , so we start by considering a quartic generating function H2 = Pi J 2 + 
foJi + PzJi-h- This does not affect $1, while <& 2 is changed into <J> 2 = <!>2 + ^ < f ) 2 
where, according to our general formula (3.15), 

(5$ 2 = {D${yP{pH 2 ) = 8(ai/3iJ 1 2 + a 2 /3 2 J 2 2 ) + 4(a 1 +a 2 )/3 3 JiJ2 • 

With the obvious notation $ 2 = b\ J 2 + 6 2 </| + b^J\Ji and with M the diagonal 
matrix M — diag(2ai, 2a 2 , ai + a 2 ), we can write in vector notation 

b = b + M(3 . 

It follows immediately from this that if ai<z 2 7^ and a\ ^ — a 2 , we can set 
b = by choosing (3 = — M~ l b: if a\ = but a 2 7^ 0, the terms 6 2 and 63 can 
be set to zero, while b\ =61, while for a\ 7^ and a 2 = we can set to zero b\ 
and 63 while 6 2 = 6 2 ; if a\ = — a 2 7^ 0, we can set to zero b\ and 6 2 , and 63 = 63. 
Finally, for a\ = a 2 = the change of coordinates has no effect on $ 2 . 
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Let us now consider this setting in the frame of Landau theory, say with a 
control parameter A £ R and with a phase transition taking place at A = (for 
definitencss, with the origin a stable fixed point for A < 0). The coefficients 
(oi, 0,2; &i, b 2 , 63) are now not constants, but in principles all depend on A; as 
customary in physical considerations, we can assume the coefficients of higher 
order terms (i.e. the 6j) do not actually depend on A. On the other hand, a 
dependence of the a, on A is physically essential. 

For a phase transition to take place, it is needed that at A = at least one 
of 01, a 2 vanish; indeed the matrix of second derivatives of $ at (0, 0) is simply 
V = diag(2ai, 202). Thus the above short discussion has an obvious physical 
meaning: we can drop a fourth order term like biJf only if the corresponding 
quadratic term diJi is positive definite throughout the range of control parame- 
ters we are considering; if the phase transition corresponds to a, passing through 
zero, we can not eliminate the corresponding b% Jf term. <0> 

We will now consider in detail example 5 (and in the next section, example 6). 
The computations are elementary but already involve some lenghty intermediate 
formulas, which are therefore not shown. 

We would like to stress that here we are only interested in discussing which 
terms can be set to zero in the Landau potential; the discussion of section 3 
also shows how one could obtain - if desired - the final coefficients as explicit 
functions of the initial ones. 

Example 5 (continued) . Recall that for Example 5 we had <f> = $1 + $2 with 
$/c homogeneous of degree 2k, given explicitely by <&i = ai J\ + aiJi + a 3 J 3 and 

$2 = b\J\ + b 2 J 2 + Mf + c x J x .J 2 + C2J1J3 + c 3 J 2 J 3 . 

The J-gradients D$/c are easily computed; these are (-D$i) = (01, 02, »3) T , and 
(L>* 2 ) = (2&iJi + ci J 2 + c 2 J 3 , ci Ji + 2b 2 J 2 + c 3 J 3 , c 2 J x + c 3 J 2 + 2b 3 J 3 ) T . The 
"P-matrix is also immediately obtained from these: 



4Ji 





2J 3 





4J 2 


2J 3 


2J 3 


2J 3 


J1 + J2 



V = 

Let us now consider an H quadratic in the J; we write this as 
H = f3iJf + (3 2 J^+ j3 3 JI + I1J1J2 + I2J1 J3 + J3J2J3 
Thus the vector (DH) is given by 

(2/?iJi +71J2 +72^3 
2P2J2+I1J1+J3J3 
2^3 +I2J1 +73^2, 

Recall that by (3.14), $1 will be unaffected, so we only have to compute 
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(5$ 2 = (D§i)V(DH 2 )- With straightforward algebra we get 

(5$ 2 - (8ai/3i+a 3 72)J 1 2 + (8a 2 /3 2 + a373)^2 2 + 
+(4(oi + a 2 )/?3 + 2a 3 (72 + 7s)] J3 + 
+ [4(ai + a 2 )7i + a 3(72 + 73)Ri^2+ 
+ [(6ai + 2a 2 )72 + 2a 3 (2/3i + /3 3 + 7i)]^i^3+ 
+[(2oi + 6a 2 )73 + 2a 3 (2/3 2 + /3 3 + 7i)]^2 J 3 • 

Hence, in the new coordinates (after the change of coordinates with generator 
H 2 ) we get 

$ 2 = $ 2 + <5$ 2 = A Ji + h Jl + .fs J 3 2 + h-h-h + T5J1J3 + h-h-h ■ 
We write 

hi fori = 1,2, 3 t f fa fori = 1,2, 3 



M 



** ^ Ci- 3 fori = 4, 5, 6 ' v ' "" \ -,; ..:•, For/ =4,5,6 
with these, and defining the matrix 

/801 a 3 \ 

8a 2 a 3 

4(ai + 02) 2a 3 2a 3 

4(ai+a 2 ) a 3 a 3 

4a 3 2a 3 2a 3 Qa x + 2a 2 

V 4a 3 2a 3 2a 3 2oi + 6a 2 7 

the transformation is written in vector form as 

7 = / + M£ . 

We should then investigate if some component of the vector / can be set to 
zero by a suitable choice of the vector £. 

This is a standard linear algebra problem, but the relevant point here is 
that this depends on the range of values assumed by the af, it depends on the 
vanishing or otherwise of the a» in the control parameter range of interest for 
the transition under study. 

In particular, the determinant of M is 

\M\ = 2 10 (01 + a 2 f (4aia 2 - a§) (3a? + 10aia 2 + 3a^ - a\) , 

and if this is nonzero we can set / to zero by choosing £ = —Mf. Note that 
here the basic invariants are not independent, which explains why it may be 
possible to eliminate all the quartic terms even when one of the aj's vanishes - 
provided the other two satisfy some nondegeneracy condition which can be read 
off the above formula for \M\. 

We recall again that quartic terms can be eliminated without affecting the 
physical validity of the Landau potential only if the latter remains thermody- 
namically stable. ^> 
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5 The Sergienko-Gufan-Urazhdin model 

In this final section we apply our discussion to the model which originated this 
work, i.e. the Sergienko-Gufan-Urazhdin model for highly piezoelectric per- 
ovskites [9] ; this is an example of a model too complex to be easily dealt with 
on the basis of "semi-intuitive" considerations for the simplification of the Lan- 
dau potential. We just consider reduction of the Landau potential according to 
our general discussion, without any consideration on the analysis of the resulting 
potential (for this the interested reader is referred to [9]). 

From the point of view of the present paper, the SGU model falls in the 
group-theoretical case considered in Example 6. We recall that in that case 
(with a slight change of notation for what concerns the coefficients appearing in 
$) the basic invariants are 

J 2 i 2 i 2 t 22i 2 2 i 22 t 2 2 2 

i = x + y + z , J2 = x y + y z + z x , J3 = x y z ; 

hence the "P-matrix is 

V = 4 




We write the most general Landau polynomial in the form 

6 

• 7 '/, 



$ = ^* A 



where the \?fc are homogeneous of degree 2k; these are, explicitely, 



*2 =hJ2 + b 2 Ji 

*3 = ci J 3 + c 2 Jf + C3J1 J 2 

#4 = di J? + d 2 Ji + d 3 J?J 2 + d 4 Ji J 3 

*s = A 4 + hJlh + h-Jl J 3 + /4J1 4 + h-hJz 

* 6 = .91 J? + giJ 2 + gs Ji + giJ?J 2 + gnJiJz + 9$J\Jl + giJiJiJ* 



Note that at the phase transition where the full G symmetry is broken, neces- 
sarily a(Ao) = 0: thus we cannot accept transformations with coefficient where 
a enters with a negative exponent. 

In [9], Sergienko, Gufan and Urazhdin argue that one can consider the re- 
duced Landau polynomial 



* = 71 -A + 72^2 + 73^3 + r x j\ + r 2 jf + r 3 J| ; 

in our present notation this means setting all coefficients to zero at the exception 
of {a, 61,62, ci,d 2 ,g 3 }- 

This statement is based on the observation that (with the notation used 
in the present note) all the other terms lie in the range of the V matrix, as 
explained in appendix B of their paper [9] . 
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Note that they do not discuss the case where some of the lowest order coef- 
ficients vanish: we have seen that this corresponds to working at a fixed value 
of the control parameter, arbitrarily near but however bounded away from the 
phase transition point. 

5.1 Reduction of the Landau polynomial with the Lie- 
Poincare algorithm 

We assume that both coefficients b\ and 62 are nonzero. 

We start with a generating function Hi G $4, written in general as 

Using our general formula (3.15), we get 

<5* 2 = {D^ 1 )V{DH 2 ) = 8(a/? 2 Ji 2 + a/3iJ 2 ) ; 

thus around the phase transition (a = 0) this vanishes and cannot be used to 
simplify the Landau polynomial. In other words, we get 

$ 2 = *2 + <5*2 = (&i+8a/3i)J 2 + (6 2 + 8a/3 2 ) Jj 2 := \ J 2 + 62J?; 

to set hi = we should choose fli = — 6j/(8a), which is singular at A = Ao since 
o(Ao) = 0. 

5.1.1 Terms of order six 

Consider the effect of this same change of coordinates on the term ^3: we have 

$ 3 := * 3 + <5*3 = (ci + 126i/3i) J 3 + 

+ (c 2 + 46i/?i + 166 2 /?i + 166i/3 2 ) Ji J 2 + (c 3 + I662A) J? := 
:= ci J 3 + c 2 Jf + c 3 Ji J 2 • 

Thus we have 

c = c + M c j3 

with the matrix M c defined as 

36i 

M c := 4 ( (fe!+46 2 ) 46i 
46 2 , 

For generic 6^., we can set to zero any two of the three Cfe. Once we have chosen 
these, and therefore we have determined the values of f3\ and /?2, it is immediate 
to compute the change also in higher order coefficients, i.e. determine also the 
(di,fi,gi). We stress that here "generic" means in particular that we need 
61 / ^ 6 2 , and h ^ -46 2 . 
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More generally, we can eliminate any two variables provided the determinant 
of the corresponding two-dimensional submatrix is nonzero: thus we have the 
following table of conditions for the different choices to be possible: 

ci = c 2 = 61 ^ 

ci = c 3 = 6ife 2 ^ 

c 2 = c 3 = (61 + 46 2 ) 6 2 7^ . 

5.1.2 Terms of order eight 

Let us now consider the effect of the change of coordinates with generator 

#3 = 7i^3 + 72^1 + 73^1^2 ■ 

We will still write the coefficients of *& as a, 6j,Cj, ....; however we stress that 
these will be different from the original ones, as they have been changed in the 
first step, i.e. under the change of coordinates generated by H 2 : thus it would 
be more precise to write them as (a, 6j, ci, di, ...), and we omit the hat for ease 
of notation. 

We have of course a 6^3 proportional to a and thus of no use near the 
transition point. As for (5^4, we get 

5*4 = (24627a) ^1 + (86173) Jl 

+ (246 l72 + 46173 + 246 273 ) J\J 2 
+ (86171 + 24& 2 7i + I26173) Ji J 3 ; 

hence \&4 = d\Jf + d 2 J 2 + ds J\ J 2 + di J\ J3 with, in vector notation, 

d = d + M d 7 
where the matrix Md is given by 



M rl 



Considerations similar to those applying for H 2 and ^3 apply also here. 

Note that, due to the structure of the matrix Md, for 61 7^ 7^ 62 and 
61 7^ 662, we can eliminate any two of d\,d 2l d 3 , and moreover (if 61 7^ — 36 2 ) 
also di. We have the following table for conditions required for the different 

choices' 

d 2 =d 3 =d 4 = bj(h+ 36 2 )^0 

d 1= d 3 =d 4 = 62(61 +66 2 )(6i +36 2 ) 7^0 

d 1 =d 2 =d 4 = 6162(61 + 362)7^0 






66 2 











26i 





661 


61 + 66 2 


(26i+66 2 ) 





36i 
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5.1.3 Terms of order ten 

We pass then to consider change of variables with generator 

H 4 = intf + % J 2 + V3J1J2 + V4J1J3 ■ 

As usual, the lowest term on which it acts is ^4, but the change is proportional 
to a and thus vanishes at the critical point. Let us then look at its effect on ^5, 
which will be mapped to 



*5 = JlJl + I2JU2 + hJlh 



/ 4 JlJ 2 



/5^2^3 



In this case we get 
5*5 = 



hence we can write 



(326 2 77l) Jf + (32&1T71 + 461773 + 326 2 %) J? Ji 
+(1261773 + 861774 + 326 2 77 4 ) J? J3 
+ (861772 + 3263772 + I661773) Ji J 2 2 
+ (2461772 + 861 774 )J 2 J 3 ; 



/ = / + M fV 



where the matrix Mf is given by 



M f 



/8& 2 
36i 



V 








(2&i + 86 2 ) 
661 



(&i 





+ 8: 

36i 

46i 




(&i 






+ 8 


26i 



/ 



h-- 


--H-- 


= .U = 


--h 


h-- 


--h-- 


= U = 


--h 


h-- 


--H-- 


= U = 


--.h 



The usual considerations apply here: we can in general eliminate four out of the 
five coefficients, and the table giving conditions for each choice is the following: 

6i(36i + 206 2 )^0 

6i6 2 (36i + 206 2 )^0 

6162(61 +46|)(6i+ 86i) 7^0 

/i = J2 = / 3 = /5 = 6162(61 + 862)^0 

h = h = h = A = 62(61 + 46|)(6i + 86D ^ 

5.1.4 Terms of order twelve 

Finally, let us consider a change of coordinates with generator 

H 5 = #! Jl + tf 2 Ji 3 J 2 + ^ 3 Jl -h + VaJiJI + V5J2J3 ; 

again we write the coefficients of the $>k with no hat for ease of notation. As 
we know, the only two terms affected by this will be $5 (with 6Q5 proportional 
to a and thus of no use) and 

*6 = .gi Ji 6 + g 2 Jl + g 3 Jl + g^l-h + gsJ^Js + g% Jl Jl + grJi J2 J 3 ; 
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for the latter we have 



(5$ 6 = (4062tfi)Jf + (86n?4)J| + (126i^ 5 )^ 2 



3 



+(40Mi + 4&102 + 406 2 ^ 2 ) Ji 4 J 2 
+(126i<? 2 + 86i<? 3 + 406 2 tf 3 ) J? J 3 
+(246ii? 2 + 861 1? 4 + 406 2 i? 4 ) J? J| 
+ (166i7? 3 + 24M 4 + 126i<? 8 + 4062^5) J1J2J3 

Thus, with by now standard notation, we get 



M g t 



where the matrix M g 



M g = 



/106 2 \ 

26i 

3&! 

IO61 (61 + IO62) 

36i (26i + 106 2 ) 

661 (26i + 106 2 ) 

V 4&i 661 (36i + 106 2 )7 

We can in general eliminate five out of the seven coefficients; due to the sparse 
nature of M g , however, we have several limitation. A complete table giving 
conditions for each allowed choice would be rather long, and is easily built by 
considering the determinant of the corresponding submatrices; it is thus omitted. 
In view of comparison with [9], it is however interesting to consider the 
choices which do not set gs to zero; for these to be possible the determinant 
of the corresponding submatrix must be zero, and hence we have the following 
table of conditions for each choice to be possible: 

.92 = .94 = .95 = .96 = .97 = 61(61 + 56 2 )(36i + 106 2 ) 7^ 

.9i = .94 = .95 = 96 = .97 = 62(61 + 56 2 )(6i + 106 2 )(36i + 106 2 ) ^ 

.91 = .92 = .95 = .96 = .97 = 6162(61 + 56 2 )(36i + 106 2 ) 7^ 

<7i = 52 = .94 = .96 = .97 = corresponding determinant is zero 

.91 = .92 = .94 = .95 = 97 = 6162(61 + 56 2 )(6i + 106 2 )(36i + 106 2 ) ^ 

gi = 52 = .94 = .95 = .96 — corresponding determinant is zero. 

Needless to say, when the corresponding determinant is zero, the choice is never 

allowed. 

5.2 Comparison with the paper by Sergienko, Gufan and 
Urazhdin [9] 

In [9] the SGU model is discussed for a fixed (noncritical) value of the control 
parameter, identified with a in our notation. It is claimed that one can always 
set to zero all coefficients at the exception of six; in the present notation these 
correspond to a, 61, 62, Ci, d,2, 53. The statement of [9] is based on the observation 
that all the other terms lie in the range of the V matrix. 
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Here we have discussed the SGU model for a in a range of values including the 
critical value a = 0. Hence we required that only terms which lie in the range 
of the matrix V when this is applied to vectors which are uniformly different 
from zero throughout a range of parameters including those at which the phase 
transition takes place, are eliminated. It is thus not surprising that we get 
results which are not coinciding with those of [9]. 

As mentioned above, the SGU approach corresponds to working at a fixed 
value of the control parameter, arbitrarily near but however bounded away from 
the phase transition point; needless to say, this can also provide a wealth of useful 
informations about the behaviour of the system near the phase transition. 

Thus, the present comparison is just aimed at showing to which extent one 
can reproduce the result of [9] by the simple and straightforward method exposed 
here. On the other hand, it also shows which of the possible extra terms should 
be checked with the more powerful methods of singularity theory, i.e. for which 
extra terms in Landau potential one has to check the qualitative behaviour is 
not affected. 

Let us now look in detail at the SGU reduced Landau potential from the 
point of view of the discussion conducted in the present section. (Recall we 
assume throughout that b\ ^ 0, b 2 ^ 0.) 

Summarizing our discussion for the SGU model with varying control param- 
eter, we have shown that if b\,b 2 do not satisfy a set of "resonance relations" 
(detailed above), then we can always reduce to a case where the only nonzero 
coefficients in the Landau polynomial, beside a, b\ and b 2 , are: one of the Ci, 
one of the d\,d 2 ,dz, one of the fi, and two of the gi. 

If &i, b 2 satisfy some of the "resonance relations" given above, then some of 
the coefficients can not be set to zero; see the above discussion for the different 
cases. Let us look more closely to these from the point of view of comparison 
with [9], always assuming b\ and b 2 are nonzero (this condition will not be 
mentioned in the following). 

We have seen in our discussion of terms of order six that one can set c 2 = 
c 3 = provided b\ 7^ —4b 2 . Note this is an extra condition, not mentioned in 
[9] . On the other hand, according to our discussion other choices would also be 
possible: in particular, one can always set c\ = c 2 = or ci = C3 = 0. 

As for order eight terms, we have seen that the choice d\ = d^ = d& = Q is 
legitimate provided b\ is neither equal to Zb 2 nor to 662. Again, on the one hand 
this restriction is not mentioned in [9], and on the other hand our discussion 
shows that other choices are also legitimate. 

Coming to terms of order ten, all of these are set to zero in [9]; our discussion 
shows that - as far as the Lie-Poincare theory is concerned - at least one of these 
must be kept different from zero if we consider a critical range of values for a; 
If the bi satisfy some resonance relation, our choice is restricted. To fix ideas, 
we set f 2 = fs = fi = fc, = 0, which is fine provided 3&i 7^ — 20b 2 . 

Finally, concerning terms of order twelve, in [9] these are all set equal to 
zero except for g%. According to our discussion, and again within the limits of 
Lie-Poincare theory, one should keep at least two of them different from zero 
(more if the bi satisfy resonance relations) . It is legitimate to choose 53 as one 



27 



of these, provided the other is neither g§ nor gj; in other words, one should 
also keep either one of gi, g2, 94, 9s as nonzero. Again to fix ideas, let us say we 
choose <7i; this is legitimate provided b\ ^ —562 and 36i 7^ — 10&2- 

Summarizing our discussion, in [9] a large number of the arbitrary coefficients 
which in principles would appear in the Landau potential are set to zero on the 
basis of singularity theory considerations; we have seen that to a large extent 
this can be justified on the basis of the much simpler theory exposed here. This 
makes use only of explicit changes of variables, with coefficients which can be 
explicitely computed by solving linear algebraic equations alone. 

Our simple considerations, on the other hand, would require other two coeffi- 
cients - i.e. other two basic monomials, one of order ten and one of order twelve 
- to be also present; moreover, on the basis of them one would suspect that the 
potential analyzed in [9] experience some special, i.e. non structurally stable, 
behaviour (even at orders not higher than eight) when one of the "resonance 
conditions" met in our discussion occur. 

It should be stressed again that these differences do not entail the claim 
that there is something wrong with the analysis of [9]: we have used a specific 
simple tool to change the coordinate expression of the general potential, while 
Sergienko, Gufan and Urazhdin allowed for elimination of other terms provided 
this does not change the qualitative predictions of the model. In other words, in 
our discussion the potential is not changed (although the coordinate expression 
changes), while they admit the potential can be changed to a different one, 
provided the two are qualitatively equivalent. 

6 Conclusions 

In the Landau theory of phase transitions, one considers an effective potential 
$ whose symmetry group G and degree d depend on the Physics of the system 
under consideration. 

As a rule, one should consider as $ the most general G-invariant polynomial 
of degree d, the latter being chosen on the basis of thcrmodinamic stability 
considerations. When such a $ turns out to be too complicate for a direct 
analysis, it is essential to be able to drop "unessential terms", i.e. to consider 
a simplified potential $ giving raise to a behaviour qualitatively equivalent to 
that generated by the general one. 

Criteria based on singularity theory and employing the spectral sequence 
technique exist and have a rigorous foundation [7, 16], but are mathematically 
sophisticated and often very difficult to apply in practice. 

Here we consider a simplifying criterion stated by Gufan [8], see the intro- 
duction. We rigorously justified a closely related reduction criterion (see section 
3) on the basis of classical Lie-Poincare theory as far as one deals with fixed 
values of the control parameter(s) in the Landau potential. 

When one considers a range of values for the control parameter(s), in partic- 
ular near a phase transition, the reduction criterion has to be slightly modified, 
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as we discussed in section 4; in particular, in order to eliminate a higher or- 
der term, certain matrices must be invertible for the full range of values of the 
control parameter (s). 

It should be stressed that in many cases one is satisfied with analyzing the 
behaviour of the Landau potential for fixed values (near the transition point) 
of the control parameter (s); in these cases the "fixed parameters" reduction 
criterion, more closely related to the Gufan simplifying criterion, has to be 
used. 

In other cases one wants to be able to "follow" the critical points of the 
Landau potential as the control parameter(s) is (are) changed; in this case one 
should use the modified reduction criterion given in section 4. 

The theory exposed here does not just provide a rigorous proof of the validity 
of the reduction criteria: it also allows to make completely explicit computa- 
tions: one passes from $ to $ by a sequence of changes of variables of a well 
defined form, and depending on a finite number of constants; the values of the 
latter can moreover be explicitely computed by solving linear algebraic equa- 
tions. 

Due to these feature, it is also straightforward to explicitly compute non- 
degeneracy conditions - ensuring the method can be actually applied to the 
problem at hand - in terms of the coefficient of low order terms in the Lan- 
dau potential <E>: indeed, the nondcgcncracy conditions correspond simply to 
the nonvanishing of the determinants of the matrices which must be inverted in 
order to implement the algorithm. 

We have considered three specific cases in detail. Two of these - i.e. examples 
4 and 5 - correspond to a two-dimensional order parameter (x,y), with the 
group G consisting, respectively, of G = {I, R x , R y , R xy } in example 4, and of 
G = {I, R xy } in example 5. Here / is the identity, R x (Ry) is the reflection in 
x (in y), and R xy the reflection in both x and y. 

In the third case studied in detail, corresponding to example 6 and the 
Sergienko-Gufan-Urazhdin model, the order parameter is three-dimensional and 
the group G, with a notation analogous to the one just used, is generated by 
{R x ,R y ,R z }, i.e. consists of G = {I , R x , R y , R z , R xy , R xz , R yz , R xyz }. 

In the final section, we studied the Sergienko-Gufan-Urazhdin model in de- 
tail; in particular, we pointed out that there are some terms which can be 
eliminated at a given (noncritical) value of the control parameter(s), but which 
should be retained - according to our method - if one wants to study a full 
range of values of the parameter(s) including critical ones. 

Finally, we note that the algorithmic procedure described here can be eas- 
ily performed by means of an algebraic manipulation language (like MAPLE, 
MATHEMATICA, MATLAB...); this also means that one can effectively tackle 
the problem of simplifying rather complicated potentials in this way. 
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